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The one-dimensional (ID) tight binding model with random nearest neighbor hopping is known 
to have a singularity of the density of states and of the localization length at the band center. We 
study numerically the effects of random long range (power-law) hopping with an ensemble averaged 



magnitude {\tij\) oc 



-J\ 



in the ID chain, while maintaining the particle-hole symmetry present 



in the nearest neighbor model. We find, in agreement with results of position space renormalization 
group techniques applied to the random XY spin chain with power-law interactions, that there is a 
change of behavior when the power-law exponent a becomes smaller than 2. 
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I. INTRODUCTION 

The one-dimensional (ID) tight-binding model of 
non-interacting electrons with random nearest neighbor 
hopping has been extensively studied since Dyson's exact 
solution to density of states of one-dimensional random 
harmonic osciUators jl], which can be mapped to the 
random nearest neighbor hopping model. Mertsching 
gave a node counting scheme to study the density 
of states in a similar model||]. Since the early 80s, 
supersymmetry methods have been used to study such 
problems with randomness]^, |^. All these studies show 
that there is a singularity in the one-electron density of 
states p{e) at the band center of the form p(e) ~ 7|T^7p- 
Concurrently, as a consequence of the Thouless theorem 
^d| , there is a logarithmic divergence of the localization 
length 5 ~ \lne\. Thus, this model exhibits a critical 
point at the band center, unlike the standard Anderson 
modelQ in one-dimension with on-site disorder, where 
all states are known to be localized!^ for any finite 
disorder. This is a consequence of the particle-hole 
symmetry that is maintained in the model despite 
randomness in the hopping when only nearest neighbor 
hopping is present. 

In this paper, we investigate in whether these fea- 
tures change if the nearest neighbor hopping model is 
generalized to include long range hopping. Since these 
properties are closely related to the particle-antiparticle 
symmetry in the Hamiltonian, we constrain the long 
range hopping to a form which preserves the particle- 
hole symmetry. (Operationally, this is done by allowing 
hopping only between sites separated by an odd number 
of lattice constants.) However, the long range hopping 
terms render analytic methods such as node counting 
theorom j|] and Thouless theorem|pd| inapplicable; 
consequently we perform numerical calculations to study 
this problem instead. 

The nearest neighbor tight binding model of non- 
interacting fermions can be mapped onto a XY spin chain 
with nearest neighbor random couplings via a Jordan- 



Wigner transformation ||i3[ . (See, e.g. McKenzieQ for 
a recent discussion of the topic). The singularity in the 
density of states of the fermion problem at zero energy 
translates to a singular magnetic susceptibility at zero 
temperature in the spin problem. Fisher |l5[| , in partic- 
ular, showed how a perturbative position space renor- 
malization group (RG) scheme using the ratio of small 
to large couplings for random spin systems employed by 
several groups in the 80s in the context of random quan- 
tum spin- 1/2 chains H, |^ becomes asymptotically exact, 
and yields the correct low temperature (energy) behavior. 
Thus the correct form of the singularity at the band cen- 
ter is obtained for the fermionic model from the asymp- 
totically exact RG. The same position space RG tech- 
nique has been used recently to study the XY spin chain 
with long range (power-law) exchange J(r) cx r^"' [ p^ . 
The behavior found there suggests that the strong disor- 
der fixed point, which determines the functional form of 
low temperature magnetic susceptibility in the spin prob- 
lem, (for nearest neighbor hopping this is related to the 
density of states singularity in the fermionic model) sur- 
vives for power law exponents a above a critical value, 
which they determined to be 2 (within their numerical 
accuracy of 5-10 per cent ). 

Using the Jordan Wigner transformation on the ID 
XY spin chain with long range couplings, while yielding 
a particle-hole symmetric fermion Hamiltonian, does not 
give a pure noninteracting fermion problem with power 
law hopping [ terms such as c|ea:p(i7r ^^^^ cJjC„)cj ap- 
pear as well ]. While the exact connection between the 
fermionic and spin models is lost when further neighbor 
hoppings are added to the Hamiltonian, we can never- 
theless expect some similarities between the two models 
- XY random spin chain and free fermions with ID ran- 
dom hopping - when the two have exchanges/hoppings 
with a power-law fall off, at least for large power-law ex- 
ponents. The result for the long range spin model also 
motivates paying close attention to our model's behavior 
when a is in the vicinity of 2. 

The outline of our paper is as follows. In the next 
section, we define the model Hamiltonian. The following 
section discusses the methods and relevant formalisms 
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used in computation. Wc then present the results for 
the density of states, localization length of eigenstates, 
and the spin-spin correlation function for the equivalent 
magnetic model, obtained from various numerical calcu- 
lations. We find that the density of states retains the sin- 
gularity at the band center, but this singularity is gradu- 
ally reduced by the long range hopping. The eigenstates 
are not well described as exponentially localized states, 
but have power-law-like tails determined by a. Finally, 
the spin-spin correlation function also shows power-law 
decay where a is not too small. 

II. THE MODEL 

The model we have studied is built on a ID bipartite 
lattice, consisting of odd and even sites along the chain. 
The Hamiltonian contains fermionic hopping terms be- 
tween the two sublattices consisting of the odd and even 
sites, respectively, while hopping within the same sublat- 
tice is forbidden. Thus we have, 

<ij> 

where i and j are summed over members of the odd 

and even sublattices respectively. This bipartite property 
preserves the particle-antiparticle symmetry and leads to 
interesting properties in the density of states and localiza- 
tion length. The hopping parameters tij are random with 
zero average, defined in terms of the distance \i — j\ and 
a Gaussian random variable, as described below. Two 
models with power law fall-off of have been stud- 

ied. Model A allows hopping between every pair formed 
by picking one site from each of the two sublattices, while 
< \tij\ > decreases with distance as a power law for each 
pair Specifically, we take 



Where Wij is a Gaussian random number with distri- 
bution 

In model B, < \tij\ > is constant, but the probabil- 
ity for a hopping term to be present decreases with the 
hopping distance as a power-law. Thus 

tij = WijM{\i-j\) 

where Wij is a Gaussian random variable as in Model 
A, and M{\i — j\) is chosen randomly to be either 1 
(with probability \i — or [with probability (1 — 

\^-J^)]■ 



We expect these two models to exhibit some simi- 
larities. We have used a = 0.2 throughout this paper. 
In both models, therefore, we have only one parameter 
a which characterizes the model. A small value of a 
produces longer range hopping. 



III. COMPUTATIONAL DETAILS 

The Hamiltonian of this model is an N x N matrix for 

a system of N sites, with elements tij = when \i — j\ 
is a multiple of 2. Because of this bipartite property, 
the Schrodinger equation can be written as two coupled 
equations when number of sites A'^ is even. Let ■tjjf'^'^ 
be the wavefunction on the sublattice consisting of the 
odd numbered sites, and tpt^^^ that on the sublattice 
consisting of even sites. If the matrix h represents the 
hopping terms from even sites to odd sites, then /i^ will 
consist of the hopping terms from odd sites to even sites. 
The Schrodinger equation can then be written as two 
equations: 

h^^odd ^ eleven _ (3) 

These equations lead to 

/l/lt^«dd = £2 ^odd 

The above equations reduce the size of the matrix of 
the NxN Hamiltonian to N/2 x N/2. Further, from the 
above equations, it is clear that the eigenvalues come in 
pairs ±e, and their eigenvectors are related by 

^even ^ ^even (g-) 
^odd ^ _^odd (7-) 

A direct consequence of the above is that the density 
of states is symmetric for each realization of the disorder. 
By exploiting the particle-antiparticle symmetry, we 
gain efficiency in our numerical calculation. 

The density of states can be obtained by directly diag- 
onalizing large matrices. In this approach, the difficulty 
is that we have to use large matrices which do not have 
significant finite size effects. Since we're interested in 
the property of the singularity at the band center, which 
by particle-hole symmetry, lies at e = 0, we need to get 
enough statistics near the band center. We identify fi- 
nite size effects by comparing the results of different sizes. 

The other approach to calculating the density of states 

is a recursive method. First, the dense matrix is cut off 
at a range large enough, so that the remaining part is 
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an acceptable approximation for the original matrix of 
arbitary size. [ For finite range hopping, this method 
gives no approximation at this stage, but for long-range 
(power law hopping) this is a different cutoff scheme, with 
a somewhat different finite size effect ] . Then the matrix 
H — e is transformed into diagonal form by a similarity 
transformation, which rotates certain columns and rows 
to eliminate off-diagonal terms. The remaining diagonal 
elements are not eigenvalues, but they retain the signa- 
ture of the matrix, i.e., a positive element implies an 
eigenvalue satisfying > e, a, negative element corre- 
sponds to an eigenvalue satisfying Ei < e. By counting 
the number of positive or negative diagonal elements we 
can obtain the integrated density of states. This process 
can be continued for arbitary length with no finite-size 
effect as in direct diagonalization, until the statistical er- 
ror of the integrated density of states is smaller than our 
requirement. For nearest neighbor model, the recursion 
equation for diagonal element is: 

tt2 

> nn — 1 

Cn = -e 7 

where i7„„_i is the only off-diagonal term in the Hamil- 
tonian, Cn is the remaining diagonal element after the 
transformation. Suppose has a distribution -F'(C) as 
n — *■ cxD, then it is given by the integral equation 

/-t-oo 
5{C + e+j)P{x)F{0dxd^ 
-oo S 

where P{x) is the distribution fmiction of Hf-^.^_^. 
Although this equation is similar to Dyson's approach 
, it has not been solved analytically due to the peculiar 
argument of the 5 function. The numerical approach 
is just to generate a sequence of C„ and calculate its 
distribution. 

To take care of the finite size effect introduced by the 
cut-off, we compare results using different cut offs. Such 
a procedure allows us to ascertain the range of energies 
for which the integrated density of states has converged: 
lower the energy, longer the cut-off required (the cutoff 
required varies roughly logarithmically with energy). 
For a fixed cutoff, the energy variation of the integrated 
density of states obtained by this method mimics 
that of the nearest neighbor i.e., finite range model - 
this sets the lower bound of energies close to the band 
center for which this method is applicable for that cutoff. 

Once we have obtained the eigenfunctions by direct 
diagonalization, we can calculate the correlation function 
of the corresponding spin model. The Jordan- Wigner 
transformation transforms the XY-spin chain in zero 
external field to a half-filled band of fermions, because 
the following identity 

Y^si^Y.'^U-^i'^ (8) 

i i 



The ground state fills the N/2 states with negative en- 
ergy which we can get by diagonalizing the Hamiltonian. 
We can further calculate the spin correlation function on 
the ground state. We compute the zz part of the spin- 
spin correlation function defined by 



C{i,j) = {SfS^) (9) 

(...) denotes the expectation value in the ground state, 
while the bar on top denotes ensemble average over the 
disorder. With = clci — ^, and i ^ j, the correlation 
function can be written as 



C{i,j)^(c]4c..c,)~- (10) 



where we have used the fact \ clcij = ^ due to half fill- 
ing. The Jordan- Wigner transformation contains a phase 
factor due to the anticommutation relation of fermion 
operators on different sites, but the phase factor coming 
from different fermion operators cancels each other in our 
expression of C(i,j). In the random singlet phaseQ ap- 
propriate for describing the nearest neighbor only model 
in the spin operator language, all these correlation func- 
tions have the same dependence on distance. The four 
fermion term can be expanded into a Hartree term and 
a Fock term. 

(cJcJc^Cj^ = V^n(»)V'rU«)V'm(j)V'm(i) 
n,m,occupied 

-^i(OV'i(j>m(j>m(z)) 

Here the sum is over all occupied states, i.e. all N/2 
states with negative energy. The Hartree term can be 
evaluated directly to be 1/4 because the wavefunctions 
are orthonormal using equations (^ and (0). Therefore 
the correlation function is given by the Fock term only. 
Further, if i and j are on the same sublattice, we see 
from equations (||) and (||) that tpnii) and ipn{j) ( n is 
among N/2 occupied states) are both eigenvectors of ei- 
ther h^h or hh"^ , and consequently orthogonal to each 
other. Thus, when i and j are on the same sublattice, 
the correlation function is exactly zero. In the numerical 
results presented below, we only display the spin cor- 
relation functions between 2 sublattices, and C{i,j) is 
replaced by C{x) where 2x + I = \i ~ j\. 



C{X) = - ^Ui)i^Ummmm{l) (11) 

m^n^occupied 

This expression can be evaluated directly, but since 
the computing time for evaluating C(x) is even more 
than that required for diagonalization, the system sizes 
we use are smaller than for diagonalization. The data 
presented for C(x) are all obtained from systems of 256 
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sites. However, we are able to discern the behavior 
reasonably well from the data for C{x) subject to this 
limitation. 



IV. DENSITY OF STATES 

Our results on model A show that the density of 
states remains singular at the band center when the long 
range hoppings are present. Figure ^ shows the density 
of states p{e) obtained by diagonalization as a function 
of energy e on a double logarithmic plot for both the 
nearest neighbor model, and for the long range model 
for different values of the power law exponent a. As 
can be seen, the singularity at the band center (e — 0) 
persists at least for large a, though its magnitude clearly 
decreases. The inset to Fig |l| compares the data for 
the nearest neighbor model along with the data for the 
lowest power law (a ~ 0.6) on a linear scale, which 
gives a clearer idea of the extent of this decrease. For 
quantitative purposes, it is better to plot the same data 
as {p{e)e)^^^ vs.lne. On such a plot, shown in fig ^, the 
nearest neighbor model is supposed to lie on a straight 
line, which it clearly does. Further, the curves show 
little deviation from the nearest neighbor model as long 
as a > 3. The deviation for smaller a are consistent 
with the singularity being gradually weakened as a 
decreases; if wc fit the data with p(e) ^ ^| ^^^^^ then 
when (7 = 2, the best fit lu is about 5. As cr becomes 
lower than 1, direct diagonalization is almost incapable 
of revealing any detail of the singularity: we only see a 
large value at the band center, and the density of states 
approaches the Wigner semicirclc[^. Nevertheless, the 
thermodynamic limit remains well defined upto a = 0.5, 
below which the bandwidth starts to increase with 
system size, i.e. we need to scale the hopping magnitude 
with system size to have a sensible thermodynamic limit. 
This critical value of a can be predicted by writing 
the model in path-integral form, using either replica 
technique or supersymmetry, and averaging over the 
random variables. Such an approach gives a four fermion 
term proportional to S|i — where i and j are site 

indices, therefore when a is less than 0.5, the sum will 
diverge. 



the integrated density of states N{e) (from to e). For 
the nearest neighbor model, the exact asymptotic form 
is given by: 



Fig y plots the inverse square root of N{e) versus In e; for 
the nearest neighbor model, the expected straight line 
behavior is seen. For power law hopping, the data shows 
measurable curvature certainly for a — 2.5. For larger a 
it is difficult to see whether the data suggest curvature, 
or simply a changing slope with decreasing a. While it 
is tempting to fit these curves with a form like 

A^(6)^|ln6r^ 

, which will lead to an singularity in the density of 
states like p(e) ~ l/(e| lne|'^+^), the data are better 
fit with several A. This suggests that corrections to 
the asymptotic form may be important for power law 
hopping. 

In figures |^,^, we show the accuracy of this recursive 
method. Figure ^ shows a direct comparison of the two 
methods. The recursive method agrees with the diagonal- 
ization results before the finite size effect sets in. Figure 
^ shows the convergence of the recursive method as the 
range of cut-off is increased. (Below a certain energy, 
which corresponds to distances beyond the cutoff length, 
the recursive method will behave like nearest model, since 
the wavefunction spreads out of the range of hopping) . 

The density of states allows us calculate the specific 
heat and spin susceptibility for this this model of non- 
interacting fermions. Using the standard formulae [p^, 
we can see that the specific heat prefactor (7, where 
Ct, = 7T) and susceptibility (x) at low temperature are 
singular at the center of the band. Thus, in the vicinity 
of the band center the zero temperature susceptibility is 
the form: 

1 

X ^ TjlnTl'^-i 

if the singularity in density of states is ^^j^Tp- A sim- 
ilar formula holds for 7. 



We have diagonalized several different sizes of samples 
to ensure that the data shown arc not corrupted by finite 
size effects. Fig. |3| shows an example of a finite size 
effect, which appears as a size-dependent rounding of 
density of states at low energy. All the data plotted cor- 
respond to N = 1024 sites unless stated otherwise, and 
do not suffer significant finite size effects. The results of 
model B are very similar to model A, except that the 
smallest a for a proper thermodynamic limit to exist is 1. 

Fig H shows our results for the density of states ob- 
tained using the recursive method. Here, one obtains 



V. LOCALIZATION LENGTH 

The nearest neighbor model is known to have a sigu- 
larity in the localization length (^) at center of the band. 
It's asymptotic form is 

f^llnel (12) 

This singularity can be deduced from the singularity 
of density of states using the Thouless theorem [pT|. 
However, in the longe range hopping model, the theorem 
does not apply, and our numerical calculation suggests 
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that the behavior of the two models is rather different. 

For finite energy, all states are found to decay from a 
central maximum in both models, and the decay becomes 
weaker as the band center is approached, as in the nearest 
neighbor model. However, the behavior with distance is 
rather different for the two models. 

• In the case of model A which has genuine long range 
(power law) behavior of the hopping parameter tij , 
the tail of the wavefunctions actually decay in a 
power-law manner, iJj{x) x~'^ . This form is obvi- 
ously determined by the power-law long range hop- 
ping term. If we apply the usual method of looking 
at the asymptotic behavior to determine localiza- 
tion length, the localization length is infinite for 
any power law exponent! 

• In the case of model B, the wavefunction is found 
to be decaying exponentially at long distances, like 
the nearest neighbor model, and the localization 
length can be obtained by several methods, which 
agree with the theoretical prediction. 

Fig shows a log-log plot of the averaged probability 
density (wavefunction amplitude squared IV'P) a-s a 
function of the distance from the center of the wavefunc- 
tion, averaged over typically 512 states, as a function of 
energy away from the band center for model A. At long 
distances, the behavior is clearly linear on this double 
logarithmic plot, implying a power law decay at long 
distances. Fitting the data shows clearly that the decay 
is related to the power law behavior of tij - the tail of 
IV'P decays as \ ^_^ly\~'i„ ■ This tells us that we cannot 
use the localization length, as defined by Thouless for 
an exponentially decaying state here. However, we may 
usefully define moments of the wavefunction to compute, 
for example, inverse participation ratios |p6[. 

The problem of defining the localization length in 
model A is not shared by model B, where we find that 
the wavefunctions always decay exponentially. To show 
this difference, in Figures || to |ll| we have plotted typical 
wavefunctions for model A and model B (on a log-linear 
plot of probability density versus distance) near the cen- 
ter and in the tail of the band. The difference between 
the two models near the center of the band is not obvious 
due to the large fluctuations, but in the tail of band, they 
look clearly different - model B shows straight exponen- 
tial decay (see Fig. down to 50 orders of magnitude 
for while model A (see Fig. ^ ) shows clear upwards 
concave curvature in the tail, characteristic of a slower 
decaying function (like a power law) at long distances. 

VI. SPIN CORRELATION FUNCTION 

We now present results for the spin-spin correlation 
function for the associated model in terms of the spin 



operators obtained by a Jordan Wigner transformation. 
The motivation is to compare the results with those ob- 
tained for the long range random antiferromagnetic XY 
spin chain[|2). We reiterate that because of long range 
hopping, our model contains terms in addition to those 
in the pure power law XY model studied by Houck and 
Bhatt; however, because of the existance of the same long 
range two-spin interactions in both studies, there may be 
several points in common. 

Figure |l^ shows the average correlation function C {x) 
obtained by averaging 256 samples of model A each hav- 
ing 256 sites plotted on a double logarithmic plot. All 
the values of C{x) are negative (corresponding to an- 
tiferromagnetic correlations). Our results suggest that 
there are three regions of distinct behavior, which can be 
summarized as given below: 

• For fast power law decays {i.e., exponents a > 2 ), 
the long distance behavior of the correlation func- 
tion remains unchanged from the nearest neighbor 
model. Thus in Figure the curves are parallel 
to each other at large x values, consistent with the 
result for the nearest neighbor model, for which a 
slope of 2 is predicted ||l5|. (For example, a best fit 
of C{x) for X within the interval [20, 60] a = 2.2 
yields a slope of 1.97). 

• When a gets below 2, the slope begins to change. 
For a close to but less than 2, the slope appears 
to be given by a itself, implying C(x) ^ x~'^ . ( A 
least squares fit to C{x) of the form C{x) — 
gives V = 1.79 for a — 1.8, v — \.AA for a = 1.5 
and V — 1.12 for cr — 1.3). Our numerically ex- 
act results thus show that the model's magnetic 
counterpart changes its low energy behavior differ- 
ent behavior for power law exponents a below 2. 
This is precisely the value around which Houck and 
Bhatt iQ found that the perturbative real-space 
RG procedure appears to break down, perhaps sig- 
nalling a change of phase. 

• Below (7 1, the deviations of C{x) from the near- 
est neighbor model become rather significant. To 
exhibit this more clearly, we show a linear plot in 
Figure C(x) seems to decay very slowly, and 
our data are consistent with it approaching a finite 
(negative) value, implying antiferromagnetic order. 
(For example, the curve for a = 0.7 is well fitted 
by + rf, where v = 0.49 is 1/2 within numerical 
accuracy). We caution, however, that for such low 
power law decays, finite size effects can be large, 
and more detailed calculations with larger system 
sizes and more samples is necessary before this re- 
sult can be stated with certainty from numerical 
studies. 

In summary, the numerical results on the spin correla- 
tion function from the long range fermion model (which 
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we can solve numerically exactly) supports the earlier ob- 
servations on the XY chain with random long range cou- 
plings using perturbative numerical RG methods that the 
random singlet phase is unstable for power law couplings 
with exponents less than 2. However, unlike the numer- 
ical RG study, which sees this breakdown of the 
RG scheme, we are able to go into the new phase, which 
appears to be characterized by continuously varying ex- 
ponents, like a critical phase. We also find evidence for a 
possible transition to long range order at still smaller a. 
It should, however, be borne in mind that these observa- 
tions, are from numerical calculations in finite systems, 
and subject to statistical errors due to finite sampling of 
the quenched random variable. 

VII. SUMMARY 

In this paper, we have presented results of a numer- 
ical study of a one dimensional lattice model of non- 
interacting fermions with random long range (power-law) 
hopping, which maintains particle-hole symmetry of the 
nearest neighbor model by allowing hopping only be- 
tween even and odd sites (i.e., no hopping allowed be- 
tween odd sites, or between even sites). We have studied 
two models - one with genuine power-law hopping, and 
another with long range hopping with a power law fall-off 
of the probability of such a hopping to be present. The 
results on density of states, localization and spin correla- 
tion function of the two models have been presented and 
analysed. 

For the density of states, we observe that the singu- 
larity at the center of the band, present for the nearest 
neighbor model, is weakened by long range hopping 
(Model A). The change is gradual, and at least for 
power-law exponents a greater than 3, is consistent with 



a change in the prcfactor of the singularity. For a less 
than 3, though, the numerical data for the desity of 
states in the range available appear to fit better with 
a somewhat different power of the logarithm of the 
energy. Beyond u = 1, the data are consistent with 
no singularity being present, until a < 0.5, when the 
thermodynamic limit becomes ill-defined. Similar results 
are seen in Model B, except that the thermodynamic 
limit becomes ill-defined at cr = 1. 

The two models exhibit rather different behavior of 
the electronic wavehmctions. Model B is conventional, 
in that its wavefunctions are exponentially localized, just 
as the eigenstates of nearest neighbor model. In model 
A, however, the wavefunctions are actually localized in 
a power-law manner rather than exponential. Conse- 
quently, the usual definition of localization length in 
terms of the logarithm of the long distance behavior of 
the wavefunction is invalid; however, several inverse par- 
ticipation ratios can still be defined. 

By transforming the fermion model back to a spin 
model using Jordan- Wigner transformation, we calcu- 
lated the spin correlation function in the ground state. 
Based on the data, three different phases may be pos- 
sible - (1) the random singlet phase, which seems to be 
stable for power law exponents down to a = 2: (2) a 
critical type phase with a continuously varying exponent 
of the power-law characterizing the spin-spin correlation 
function between a = 2 and cr = 1; and (3) a possibly 
long-range ordered phase for a < 1. 
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FIG. 1: Density of states of model A for different values of the 
power law exponent a, shown on a double logarithmic plot, 
showing the singularity at the band center (e = 0). The inset 
shows the data on a linear scale. The singularity becomes 
very weak when a is less than 2, and is not discernable for 
the sizes studied when a is below 1. 




FIG. 2: The same data as in figure |l|, plotted in a differ- 
ent way as motivated in the text. As expected, the density 
of states of the nearest neighbor model asymptotically falls 
on a straight line. The upwards bending of the curves with 
decreasing a, evident for ct < 3 shows the weakening of the 
density of states singularity by the long range hopping. 
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FIG. 3: Finite size effects on density of states, which appear 
at the low energy end of these curves. Data are density of 
states for four different size systems with the same a = 2.5, 
with length varying from N = 128 to 1024. 



by Recursive Method 
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FIG. 4: Integrated density of states, computed by the re- 
cursive method, with cut-off at 100th neighbor. Within the 
displayed energy range, these curves are checked to be free of 
effects due to the finite cut-off in the hopping range. Note 
that the curves are smoother than the results obtained from 
straight diagonalization. The nearest neighbor model is ex- 
pected to asymptotically fall on a straight line with a slope 
equal to -1. Curves fitted by the formula N(e) ~ l/jlnel", 
yield v somewhat larger than the value 2 appropriate for the 
nearest neighbor model, for a below 3. 
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FIG. 5: Comparison between the recursive method and the 
direct diagonalization method for the density of states. The 
solid curve is by direct diagonalization, while the points are 
density of states obtained by taking the difference of inte- 
grated density of states from the recursive method. The agree- 
ment between the curves is good, and no systematic errors are 
found. 
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FIG. 6: Convergence of the recursive method exhibited using 

different cutoff jjaramctcrs for the hopping. Small values of 
the cutoff lead to deviations starting from higher energies. 
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FIG. 7: Double logarithmic plot of the probability density in 
real space for individual eigenstates plotted from the center of 
each state, averaged over eigenstates with a given energy e for 
model A with a = 3.5. Small values of e correspond to states 
near the band center. We can see as the energy approches 
the band center, the states become more dclocalizod. In the 
tail, for all energies, the profile of the wavcfunction actually 
decreases in a power-law fashion, with an exponent exactly 
determined by the value of a. Similar plots for other values 
of a in model A show the same feature. The inset with the 
probability density on a semilog plot, shows distinct curvature 
in the tail, unlike an exponentially localized state. 
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FIG. 8: A typical wavefunction of model A near the band 
center(small e) shown on a log-linear plot 
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FIG. 9: A typical wave function of model A away from the 
band center(larger e) on a log-linear plot. The tail shows 
clear curvature on this plot, implying a functional form that 
decays slower than exponential. 
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FIG. 10: A typical wavefunction of model B away from the 
band center. A clear exponential decay of the wavefunction 
amplitude over 50 decades is seen, as in the nearest neighbor 
model. 
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FIG. 11: A typical wavcfunction of model B near the band 
center. Despite strong fluctuations, the main profile is basi- 
cally exponential, as in the nearest neighbor model. 
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FIG. 12: Spin-spin correlation function in the ground state 
(T = 0), for the corresponding spin model obtained using a 

Jordan- Wigner transformation. The bold straight line, which 
decays as is a gtiidc for the eye. Wo sec that when a is 
above 2, including nearest neighbor model, C(.t) exhibits this 
inverse square behavior. Below a = 2, C{x) is better fitted by 
More sample averaging is necessary to smooth the noisy 
tail at long distances. 
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FIG. 13: Spin-spin correlation function, as in figure |l^, plot- 
ted on a linear scale. For a < 1 C{x) decays very slowly, and 
seems to approach a constant value at large distances. 



